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Abstract 

The ANTARES high energy neutrino telescope is a three-dimensional array of about 900 photomultipliers 
distributed over 12 mooring lines installed in the Mediterranean Sea. Between February and November 2007 
it acquired data in a 5-line configuration. The zenith angular distribution of the atmospheric muon flux 
and the associated depth-intensity relation are measured and compared with previous measurements and 
Monte Carlo expectations. An evaluation of the systematic effects due to uncertainties on environmental 
and detector parameters is presented. 

Keywords: atmospheric muons, neutrino telescope, depth-intensity relation 



1. Introduction 

The ANTARES collaboration has recently completed the construction of a neutrino telescope in the 
Northern Hemisphere. Its location is about 40 km off the coast of Toulon, France, at a depth of 2475 m 
in the Mediterranean Sea (42°48'N, 6°10'E). The main goal of the experiment is the search for high-energy 
neutrinos from astrophysical sources such as Active Galactic Nuclei, microquasars, supernova remnants and 
gamma-ray bursters. A neutrino telescope in the Northern hemisphere includes the Galactic Centre in its 
field of view and is complementary to the AMANDA and IceCube Antarctic telescopes [l| . Neutrinos are ob- 
served indirectly via the detection of the Cherenkov light emitted along the trajectory of high energy charged 
particles emerging from neutrino interactions. The Cherenkov light is collected by light sensors (photomul- 
tipliers - PMTs) distributed in a three-dimensional array. The PMTs are arranged on twelve vertical lines. 
The first detection line was deployed in February 2006 and a second became operational in October of the 
same year. Three more lines were connected in January 2007 and five more in December 2007. The ap- 
paratus reached its complete configuration when the last two lines were deployed and connected in May 2008. 

The telescope is optimised for the detection of muon neutrinos, since muons resulting from charged 
current interactions can travel kilometres and are almost collinear with the parent neutrinos at energies 
above 1 TeV. The technique relies on discriminating upward-going neutrino-induced muons from the flux 
of downward-going atmospheric muons, which represent the majority of events in the ANTARES detector. 



* Corresponding author 
Email address: Aimarita.Margiotta@bo.infn.it (A. Margiotta) 
Also at University of Leiden, the Netherlands 
On leave at DESY, Platancnallcc 6, D-1573S Zcuthcn, Germany 



Preprint submitted to Astroparticle Physics June 28, 2010 



Atmospheric muons are produced mainly by the decay of charged pions and kaons resulting from the inter- 
action of high energy cosmic rays with atomic nuclei in the atmosphere. Although they are a background 
for neutrino detection, the atmospheric muons are useful to verify the detector response. 

This paper reports on the determination of the zenith angular distribution and on the measurement of 
the flux of atmospheric muons using data taken between February and November 2007 when the detector 
comprised 5 lines. The experimental results are compared with atmospheric muon flux simulations, with 
previous ANTARES measurements obtained with the first detector line [2J and with the results of a low 
energy event selection [3j] . 

2. Experimental setup and data selection 

The full ANTARES detector consists of 12 flexible lines, each with a total height of 450 m, separated by 
distances of 60-70 m from each other. They are anchored to the sea bed and kept near vertical by buoys at 
the top of the line. Each line carries a total of 75 10" Hamamatsu PMTs housed in glass spheres, the optical 
modules (OM) [J], arranged in 25 storeys (3 optical modules per storey) separated by 14.5 m, starting 100 m 
above the sea floor. Each PMT is oriented 45° downward with respect to the vertical. A titanium cylinder 
in each storey houses the electronics for readout and control, together with compasses and tiltmeters used 
to measure the heading and the inclination of the storeys. To reconstruct the neutrino direction with high 
precision a good knowledge of PMT position and of hit arrival times is necessary. The positions of the 
optical modules are measured by a system of acoustic transponders and receivers distributed over the lines 
and on the sea bed. A system of LED beacons housed in some of the storeys and a laser beacon located at 
the bottom of two of the lines are used for timing calibration |5( . The PMT signals are digitized by a custom 
built ASIC chip. For analog pulses which are larger than a preset threshold, typically 1/3 photoelectron 
(p.e.), the arrival time and the integrated charge of the pulse are measured (referred to as L0 hit) and the 
corresponding data are sent to shore. The data stream is processed by a computer farm in the shore station 
which searches for differentphysics signals according to predefined trigger conditions. The DAQ system is 
described in detail in Ref. [6j. 

During the 5-line data taking period, the trigger algorithm required at least 5 causally connected LI 
hits, where an LI hit is defined either as 2 L0 hits in coincidence within 20 ns in two optical modules on the 
same storey or as a single L0 hit with an amplitude larger than 3 p.e. 

In Figure [T] the typical counting rate registered by an optical module is shown. It is characterized by a 
minimum constant rate, the so called baseline rate, originating from 4n K decay and bioluminescent bacte- 
ria. Superimposed on the baseline are occasional bursts, of a few second duration, attributed to luminous 
emission by macro-organisms. 

For this analysis, only runs fulfilling the following quality criteria are selected: i) at least 300 out of 375 
installed optical modules are active during a run (80%); ii) the baseline rate is below 120 kHz; iii) the burst 
fraction is less than 20%. The burst fraction is defined as the fraction of time during which the instantaneous 
optical background rate exceeds the baseline rate by at least 20%. With this selection the active detector 
time was about 90 effective days, during which more than 10 7 atmospheric muon triggers were collected. 

3. Analysis chain 

The analysis chain comprises the simulation of the underwater muon flux induced by the interaction 
of the primary cosmic rays in the atmosphere, the simulation of the Cherenkov light in the instrumented 
volume and the reconstruction of the muon tracks starting from the detected or simulated signals. 



X 
<u 


700 
600 
500 




400 




300 




200 




100 



- 








■ OM1 


z_ 






■ 


OM2 


- 






■ ■ 
• 


■ OM3 










I 




~ 








~ 






:■ i 




— 






V 




— 






'■• 










% ■ 




— 






: ■ 




— 






■ • : 

■ 






i 


- - kJi* 






- , , i , , , i , , , i , , , i , , , i , , , i , , , 



11m20 



11m40 



12m00 



12m20 



12m40 13m00 

Time since Runstart (min) 



Figure 1: Typical counting rates in three optical modules of a detector line storey. 



3.1. Simulation of the underwater muon flux 

Two different methods are used to describe the atmospheric muon flux reaching the detector. 

The first approach is a Monte Carlo simulation based on CORSIKA version 6.2 [7|. Two different models 
have been considered for the description of hadronic interactions: QGSJET.01 [8| and SIBYLL [9J. Showers 
were simulated over an energy range from 1 to 10 5 TeV/nucleon and zenith angles between 0° and 85°. 
For zenith angles larger than 70° the curvature of the atmosphere is taken into account. Five mass groups 
(p, He, CNO, Mg-Si and Fe) are generated with an energy spectrum proportional to E~ 2 . More than 10 10 
air showers induced by primary cosmic ray nuclei were simulated with QGSJET.01. This corresponds to 
a detector live time ranging from few days to few years depending on the primary mass, its angle and 
energy. A smaller sample was simulated with SIBYLL, mainly to evaluate the effect of a different choice 
of the hadronic interaction model. Muons reaching the sea level with energies larger than 500 GeV arc 
propagated through sea water to the detector with the MUSIC package [10j. The Monte Carlo events are 
weighted according to two different primary cosmic ray s pec tra : the Horandel/poly-gonato model (with 



rigidity dependent knee energies) ll| and the NSU model [12|,|l3(. Both are based on experimental results, 



combining direct measurements and results from extensive air shower arrays. Note that in the energy range 
relevant for us (10 - few hundred TeV) the two all-particle spectra differ by about 30%. 

The second approach uses a Monte Carlo event generator based on a parametrized description of the 
underwater muon flux (MUPAGE) [l4|. The flux, the angular distribution, the multiplicity and the energy 
spectrum of muons are obtained from parametric formulae [15] . A number of showers corresponding to one 
month of live time has been generated. 



3.2. Light production and detector simulation 

All muons reaching the detector depth with an energy above 20 GeV are transported through the detector 
active volume using MUSIC. As the number of photons emitted along the muon trajectory is large, the 
simulation of the light production is very time consuming. For this reason, a set of tables storing amplitude 
and time of Cherenkov photons hitting the PMTs is created. They are constructed taking into account the 
water properties at the ANTARES site (light absorption length and scattering model), the characteristics 
of the optical modules (geometry and efficiency) and the position, distance and orientation of an optical 
module with respect to a given muon track or electromagnetic shower. The optical background is taken from 
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Figure 2: A downgoing track, reconstructed using all 5 lines. In each display, the curve shows the reconstructed muon track, 
i.e. the z position of the Cherenkov cone, at the x-y position of the corresponding line, as a function of time. Horizontal axis: 
time in nanoseconds; vertical axis: height above the sea-bed in metres. Crosses are hits in a time window of 3 microseconds 
around the trigger; full circles are hits passing the trigger condition (see Section 2); open boxes are hits used in the final fit. 



the actual counting rate observed in the data. Dead channels and the efficiencies of the front-end electronics 
are also taken into account. 

3.3. Track Reconstruction 

A fit algorithm based on a chisquare-like minimization is used to reconstruct the tracks. It approximates 
each storey as a space point on a vertical detector line. The hits in the event are time ordered and merged 
on each storey if closer than 20 ns; in this case the time of the first hit on the storey is taken and the 
charges summed. In order to augment the weight of coincidences with respect to single high-charge pulses, 
the summed charge is further increased by 1 if the hits originate from different optical modules of the same 
storey. All hits, merged or single, having a minimal charge of 2.5 p.e. are defined as step-1 hits. Next, a 
search is made for clusters of step-1 hits by requiring the presence of two step-1 hits within 80 ns on two 
adjacent floors or within 160 ns on two next-to-adjacent floors. For lines having one or more such clusters, 
single L0 hits within the line are selected if they are causally connected to the cluster (s) of the line. The 
times and the positions of all selected hits on all lines are then used in a prefit of the track. The result of 
this prefit provides the starting values for the final chisquare minimization. This final fit yields two angles 
defining the track direction and the three coordinates of the track at a chosen time. 



Monte Carlo studies indicate that this reconstruction strategy yields a resolution on the zenith angle of 
0.7° for atmospheric muons. The assumption of a vertical detector line introduces a further uncertainty on 
the reconstructed zenith angle. The sea currents registered during the selected data period induce displace- 
ments of the OMs from the nominal position which are smaller than one metre over 350 m height. This 
corresponds to an additional error of about 0.15°. 

An example of a downward going track, reconstructed in the data, is shown in Figure [21 The horizontal 
axis indicates hit time, while the vertical axis shows the height of the fired storey in metres. On this 
picture, crosses represent hits in a time window of 3 microseconds around the trigger, full circles indicate 
hits participating in the trigger and open boxes designate those hits which have been used in the fit. 
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Figure 3: Zenith angle distribution of reconstructed tracks. Black points are the data. The dotted line is the MUPAGE Monte 
Carlo expectation. The solid line corresponds to the simulation with CORSIKA + QGS JET and the NSU model of the primary 
cosmic ray flux; the dashed line is CORSIKA + QGSJET with the "poly-gonato" model. The dashed-dotted line refers to 
CORSIKA + SIBYLL and the NSU model. The shadowed band represents the systematic uncertainty, with respect to the 
solid line. See text for details. 



4. Zenith distribution 

Figure [3] shows the zenith angle distribution of the reconstructed tracks. Each track represents a muon 
bundle, regardless of the real number of muons in the bundle. Only tracks reconstructed using hits on at 
least 6 floors are considered in this plot. The systematic errors on the Monte Carlo expectations are due 
to the uncertainties on the description of the detector, on the knowledge of the environmental parameters 
and on the primary cosmic ray composition and hadronic interaction models. The effect of the first two 
terms is shown as a grey band in Figure [3] and has been evaluated by repeating the Monte Carlo simulation 
described in Section 3.2 for various assumed input parameters. 



The values of the ANTARES environmental parameters (light absorption and scattering) have been mea- 



sured during several sea campaigns, [16I.I171 Il8l|. The uncertainty on the light absorption length in water is 



assumed to be ±10% over the whole wavelength spectrum and yields a variation of ±20% on the number of 
expected muon events. 



The detection efficiency of an optical module is a function of the angle between its axis and the Cherenkov 
photon direction. The acceptance for photons impacting the PMT at small angles with respect to the PMT 
axis (i.e. head-on) is determined with an uncertainty of ±15%. This is the region of the angular acceptance 
mostly affecting the detection of upgoing neutrino induced muons. At large angles (i.e. glancing), where 
the OMs are mostly sensitive to downgoing atmospheric muons, the uncertainty is about ±30%. The global 
effect is a systematic uncertainty of ^np on the reconstructed track rate expectation. The assumed PMT 
effective area is that provided in the manufacturer data sheets [19| . It takes into account the photocathode 
area, the collection and the quantum efficiencies. A study of 40 K decays observed in the ANTARES detector 
indicates an uncertainty of ±10% on the PMT effective area [3J and results in a variation of about ±20% 
on the expected event number. Summing the contributions from the various environmental and detector 
related systematic uncertainties leads to an overall uncertainty on the total number of reconstructed tracks 



of +45% 

01 -40%- 

The reconstructed track rate is 1.52 Hz in the data. The expected track rate with the MUPAGE event 
generator is 1.55 Hz [14j. In the CORSIKA Monte Carlo samples, the reconstructed event rates are 1.1 Hz 
for QGSJET + NSU model 0, 0.76 Hz for the QGSJET + "poly-gonato" model [Hj], 1.2 Hz for SIBYLL 
+ NSU model. There is an important dependence of the Monte Carlo prediction on the hadronic interaction 
model used to describe the development of the air showers and on the energy spectra of primary cosmic 
rays. Within systematic uncertainties the measurements cannot distinguish between the models considered. 

5. Depth-intensity relation for atmospheric muons 

The 5-line data were used to determine the depth-intensity relation of atmospheric muons, i.e. the 



vertical muon flux versus depth, h, [20j, |2 lj j . Given the muon intensity, 1(9, h ), at a certain depth h as a 
function of direction 9, the vertical intensity of the muon flux as a function of depth underwater 1(9 = 0, h) 
is calculated, using the relation: 

1(9 = Q,h) = 1(9, K) ■ \cos(9)\ ■ c corr (9) (1) 

where h — h /cos(9) is the effective slant de pth and c corr (9) is a geometrical correction factor which takes 
into account the curvature of the Earth, [22J, [23J . 

The underwater muon flux cannot be measured directly with the ANTARES detector (and in general 
with neutrino telescopes). It is calculated starting from the angular distribution of reconstructed tracks, 
multiplied by the average multiplicity of events, estimated with Monte Carlo methods, according to the 
following equation: 

N(9,h )^(9,h ) 

mh ° ) -A.„w).T.*nw (2) 

where : 

• N(9, h ) is the number of tracks reconstructed in the angular bin around cos(9). It is obtained from the 
zenith angular distribution of all reconstructed tracks and is corrected to account for the trigger and 
reconstruction efficiencies, calculated with Monte Carlo simulations. An iterative unfolding procedure 
based on Bayes theorem was applied [24j; the resulting correction is small because of the good zenith 
angle resolution compared with the bin width. 

• /j,(9, h ) is the mean muon multiplicity in the bundles, estimated from the MUPAGE simulation [15[. 
It ranges from 1.3 for vertical events to 1.1 for almost horizontal events. The CORSIKA Monte Carlo 
gives 7% smaller values. 

• A e ff (9) is the Monte Carlo computed effective area for muons at the 9 angle. 

• T and Afl(9) are the live time for the selected data and the solid angle for the cos(9) bin, respectively. 

Figure 4 shows the atmospheric muon flux, 1(9, h = 2000), computed at 2000 m depth. It is compared 
with the MUPAGE simulation (solid line). 

The measured muon flux, 1(9, h ), was finally transformed into the vertical muon flux using eq. (1). 
The results are shown in Figure [5] (black points). Statistical errors are small and not visible on the plot. 
The error band represents the systematic uncertainties discussed in Section 4. Note that other than in the 
preliminary analysis described in ref. [20j] , in this case no selection on the quality of the track fit is applied. 
Previous results published by ANTARES and a calculation obtained with a parametrization taken from 13] 
are also shown. All measurements agree within the systematic uncertainties with the model predictions. 
Several measurements from other underwater detectors are also shown for comparison. 




Figure 4: The flux of atmospheric muons at 2000 m depth (black points) as a function of the cosine of the zenith angle. The 
solid line is the expectation from the MUPAGE simulation. The grey band represents the systematic uncertainties on the 
measured flux. 
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Figure 5: Vertical muon flux of atmospheric muons for the 5 line ANTARES data (black points) as a function of the slant depth. 
Downward triangles show the results from line 1 data 01 • Full squares show the results obtained with a new method and a low 
muon energy event selection [3|. Expectation from Bugaev parametrization (dotted line) is superimposed 13]. A compilation 
of results obtained with other underwater detectors is shown: AMANDA |2j, AMANDA-II |2J| Baikal |27l. DUMAND [28JI . 
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6. Conclusions 

Using the data collected with the first 5 lines of the ANTARES detector, the zenith angle distribution of 
atmospheric muon bundles and the vertical muon flux as a function of the water depth were measured. The 
results are compared with Monte Carlo expectations obtained with a complete Monte Carlo simulation chain 
and with a parameterized calculation of the underwater muon flux at the detector depth. The major sources 
of systematic uncertainties in Monte Carlo simulations are the description of the angular dependence of the 
optical module efficiency and the theoretical models used to represent the primary cosmic ray composition 
and the hadronic interactions during the air shower development. The results of the present analysis are in 
agreement, within the systematic uncertainties, with the theoretical predictions and previous measurements. 
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